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Abstract 

The spectra of a particular class of VT symmetric eigenvalue problems has previously been 
studied, and found to have an extremely rich structure. In this paper we present an explanation 
for these spectral properties in terms of quantisation conditions obtained from the complex WKB 
method. In particular, we consider the relation of the quantisation conditions to the reality and 
positivity properties of the eigenvalues. The methods are also used to examine further the 
pattern of eigenvalue degeneracies observed by Dorey et al. in [1,2]. 
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1 Introduction 

The one-dimensional, time-independent Schrodinger equation 



with G L^{C) {i.e. is square integrable over some given contour C in the complex 

plane), is an eigenvalue problem for the energy E. Eigenvalue problems of this type are a natu- 
ral extension of eigenvalue problems on the real line [3] , and indeed, the usual requirement that 
\E'(x) G L'^ijVj just corresponds to a particular choice of contour. 

In recent years, there has been considerable interest in the study of such eigenvalue problems 
associated with Hamiltonians that have the property of VT symmetry [4-6]. The VT operator 
is an anti-linear operator corresponding to parity reflection and time reversal: 

V^{x) = $(-a;) 
T^x) = ^*{x) 



A reason for the interest in VT symmetry is that the spectrum of VT symmetric Hamiltonians 
sometimes appears to be purely real. In general, it can be shown that the eigenvalues of these 
Hamiltonians will be real, or occur in complex conjugate pairs. A connection has also been 
recognised between these types of differential equations and the field of integrable models (the 
'ODE/IM correspondence') [7-12]. 



In this paper we will consider a particular VT symmetric eigenvalue problem 



dx'^ 



6 2 ^(^ + 1) 

+ ax^ + ^ „ ^ 



^(x) = A^(a;), ^(x) G L\C) 
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This is the M = 3 case of the family of equations 
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- Uxf'' - aUx)''-' + 
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^(x) = A^(x), ^(x) G L\C) 



discussed in [1,2], with the contour C joining |x| = cxo in the Stokes sectors iS_i and Si (Figure 



]m(xj j 




Figure 1: The quantisation contour, C (M = 3) The Stokes structure is 
shown in the limit of large \x\, without showing the detailed structure near the 
turning points. 



It was shown in [1,2] that this equation, with these boundary conditions, has a spectrum which 
is: 

real if a < M + 1 + 1 2/ + 1 1 [Regions B,C,D] 

positive ifa<M+l — |2/ + 1| [Region D] . 

However, for q;>M+1 + |2/ + 1|, (Region ^4), complex energy levels may be found. For the 
M = 3 case, these four regions in the {a, I) plane are shown in Figure [21 separated by the dashed 
lines. 
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Figure 2: Regions of the a - I plane. In [1], Dorey et al proved the reality 
of the spectrum for (a, /) G BUCU D, and positivity for {a, I) & D. The dotted 
lines show a± G Z. 
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Also shown in [1,2] was an interesting structure to the boundary of the region, within A, where 
eigenvalues become complex (the 'domain of unreality'), notably the appearance of 'cusps' on 
this boundary (Figure [3]). 




0- 1 /' \ 

-8 -6 -4 -2 2 4 6 

Figure 3: Domain of unreality. The outline of the region where complex 
eigenvalues first appear, as shown in [1,2]. 

The remainder of this paper will be structured as follows: Section [2] deals with the complex WKB 
method, and how quantisation conditions derived using this method are able to account for the 
results found in [1,2]. In particular, in l2.H we give a brief summary of the complex WKB method 
and its use. Section fI72\ explains how this method can be used to obtain quantisation conditions 
for the energy eigenvalues, and in 12.31 we describe how different quantisation conditions are ap- 
plicable in this problem for different values of the parameters E, a and /. Section [2.3.11 gives an 
example of calculating a quantisation condition for some particular values of the parameters. In 
section 12.41 we demonstrate that the WKB quantisation condition approach is able to explain 
the positivity of the spectrum in the region found in [1,2]. Section [275] explains the appearance 
of complex eigenvalues in terms of the quantisation conditions, in the region where Dorey et al. 
found they could occur, [1,2] , and we conjecture how complex WKB methods might demon- 
strate reality. In section 12.61 we use the WKB quantisation conditions to calculate the positions 
of degenerate eigenvalues (and reproduce the boundary of the domain of unreality), and explain 
the formation of the cusps. Our results are compared to those found in [1,2]. Finally, Section [3] 
contains conclusions and a discussion of possible future work. 

As a brief aside, in [2], Dorey et al. defined a new set of coordinates on the {a, I) plane by 

a± = \a - M - 1 ± (2/ + 1)] 

2M + 2^ ^ 
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which we will also occasionally adopt in order to allow an easier comparison of results with [2]. 
In particular, for the M=3 case dealt with here, 

a± = ^[a-4± (2/ + 1)] 

o 

so that 

OL = 4(1 + + a_) and I = -[4(a+ — «_) — 1] 
The dotted lines shown in Figures [2] and [3] indicate a± G Z. 



2 Complex WKB 
2.1 General Principles 

For the one- dimensional, time-independent Schrodinger equation 

-—\^ + V{x)^{x) = E^{x), 

the first-order WKB approximation (see e.g. [13]) states that two solutions to this equation are 
given asymptotically for |x| ^ oo by 

'^±{x) ^exp ( ±z / {P{t))^dt] , P{x) = E- V{x). (2) 

The complex WKB method (see e.g. [14-16]) involves treating x as a complex variable, and cal- 
culating how the asymptotic form of the solutions change as they are traced around the complex 
plane. The lower bound of the integral Xq will be a zero of P{x) {i.e. a turning point of the 
classical problem), and the solutions are valid in some region of the complex plane around Xq. 
The integral in the exponent is, in general, complex valued. Because it will typically have an 
imaginary part, the solutions will contain real exponentials. If a solution contains a growing 
exponential term, it is said to be dominant, while a solution with a decaying exponential is 
called subdominant. However, if this integral is purely real, then both WKB solutions will be 
purely oscillatory, with neither dominating the other. Lines can be drawn in the complex plane, 

emanating from the turning points (zeros), marking the curves where 53m J^^{P{t))^dt) = 0. 

These are known as anti-Stokes lines. Anti-Stokes lines mark the borders between sectors of dom- 
inant /sub dominant behaviour; in a given sector one solution will be dominant, but on crossing 
an anti-Stokes line to the next sector this behaviour reverses and the solution will be subdomi- 
nant in the new sector. 

We can also find the regions where JJe \j^{P{t))2dt)\ = 0. These are known as Stokes lines, 

and are the regions where the solutions are most dominant /sub dominant. Knowing the position 
of these Stokes and anti-Stokes lines is crucial in applying complex WKB methods. 
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Beginning with a linear combination of the two solutions at one point of the complex plane, 
a globally defined solution may be obtained by following several well-established rules. The 
most important of these for our purposes can be summarised: 

1. When crossing an anti-Stokes line, the solutions 'exchange dominance', i.e. a dominant 
solution becomes subdominant, and vice versa. 

2. Upon crossing a Stokes line, the coefficient of the subdominant solution changes by an 
amount proportional to the coefficient of the dominant term. 

i.e. new sub. coefRcient = old sub. coefRcient + T x old dom. coefficient 

T is called a Stokes multiplier, and depends on the nature of the turning point that the 
Stokes line originates from. For a Stokes line emanating from a simple turning point the 
value of the Stokes multiplier T is i, for a Stokes line from a double zero, T = V2i (for 
the first-order approximation). The coefficient of the dominant term remains unchanged 
when a Stokes line is crossed. 

3. The rules given in 1. and 2. refer to a WKB solution defined in terms of a particular 
turning point crossing a Stokes/anti-Stokes line emanating from that turning point. If it 
is intended to continue a solution across a line from a different turning point, the WKB 
solution must first be rewritten in terms of this new zero. To connect solutions defined in 
terms of different turning points, Xi and X2, use: 

exp [±^ j:^{P{t))Ut) = exp (±z /-(P(t))icit) exp [±^ jl{P{t))Ut) . 

Note that it is chosen to introduce the discontinuous change in the subdominant coefficient at the 
Stokes line, so that this discontinuity is small compared to the error in the WKB approximation. 



2.2 WKB Quantisation conditions 

We study the spectrum of this eigenvalue problem via quantisation conditions obtained using 
the complex WKB method. Briefiy, this is done by starting with a WKB solution that is purely 
subdominant in the Stokes sector iS_i (Figured]). This solution can then be traced around the 
complex plane, in an anti-clockwise direction, with new subdominant exponentials appearing as 
Stokes lines are crossed. The new subdominant contributions will have a coefficient of i (the 
Stokes multiplier associated with line from a simple turning point) or \^i (the Stokes multiplier 
associated with line from a second order zero) multiplied by the coefficient of the dominant 
exponential at that Stokes line (in accordance with the established complex WKB rules - see for 
instance [14,15]). This leads to a wave function in the Stokes sector Si that has both a subdomi- 
nant and dominant component. Requiring that the coefficient of the dominant solution vanishes 
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(to fulfil the boundary conditions) leads to a quantisation condition for the energy. This is an 
approach used in [17] to investigate the appearance of complex conjugate eigenvalues in another 
VT symmetric problem. An explicit example of this type of calculation is shown in Section !^. 3. II 



2.3 WKB Quantisation conditions for this eigenvalue problem 

We now specialise our discussion to the eigenvalue problem of eq. ([I]). Complications arise in this 
problem because many different arrangements of turning points, and hence different topologies 
of Stokes and anti-Stokes lines, are found for different values of the parameters in the equation. 
This leads to different quantisation conditions being required for different ranges of values of a, 
/ and E. These different possible arrangements are summarised in Table [T] and Table [21 Table 
[T] refers to positive values of a, and Table [2] to negative. Each table is split into five rows, by 
different ranges of values of /, relative to the key values / = and / = V . The remainder of 
this section will now be used to explain the origin of this and the way the different Stokes 
structures are organised according to the parameter values. 

As energy is increased (or decreased) from 0, pairs of turning points move together, until an 
energy value is reached at which a double zero is found. The double zero then splits, and the 
two turning points move off again, perpendicular to their original path. (See Figure H]). 
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E=0 0<E<E' E=E' E'<E<E" 



Energy 



Figure 4: Changes in turning point configuration occurring with in- 
creasing energy. Example shown is for a = 3,/ = 0.5. This corresponds to 
the first row in Table [H E' is the energy value at which the zeros coalesce to 
form a double zero. E" is the energy value at which the outer of the zeros in 
the horizontal line 'line up' vertically with the outer zeros. Both E' and E" 
can be expressed as a function of a and I. For these values of a and I we have 
E' ^ 3.1075, ^ 3.5905. 
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Different Stolces structures are obtained according to where this coalescence occurs in respect 
to the other turning points, i. e. new topologies are obtained depending on whether the double 
zeros occur within, outside of, or directly in line with the other turning points. This can be 
easily seen by looking at Figure [5] (which corresponds to looking down the E = E' column of 
Table [1]); In the first entry the double zeros appear inside the other four zeros, in the second 
they line up vertically with the other zeros and in the third the double zeros are outside the box 
formed by the remaining zeros. 



O 

• • • * • • 

/</' 1=1' />/' 

> 



Figure 5: Relative positions of double zeros, a > 0, £" = E' . This corre- 
sponds to the E = E' column in TableUi 



For a given value of the parameter a, the value of / required for the zeros that coalesce to line 
up with the others, (which we will denote /'), is given by 



1 , 



/' = — + L (3) 

If / is less than this /', turning points move together within the remaining turning points (in a 
horizontal sense for positive energies and vertical sense for negative energies). For I > /', the 
turning points that coalesce are positioned outside of the others (again in a horizontal sense for 
E > and vertical for E < 0). 

Another important value of / where the behaviour changes is / = 0. For I = and / < 0, different 
sequences of arrangements occur, and these are shown in the second and third blocks of Tables 
[Hand [21 Obviously for / = (and / = —1), the order x"^ term vanishes from the potential, and 
there are now only six zeros, instead of eight. 
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Table 1: Arrangements of Turning Points, a > 
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(Figure [6]l 
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(Figure |6ll 
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E = 


< E < E'" 
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(Figure [6} 




E < ~E" 


E = ~E" 


~E" < E <0 


E = 


< E < E'r 


E = E" 


E> E" * 


Boundary of //J 


I = 

















t-£;' = -E" for I = I' 
»E' = E" for I = I' 

* Change in Stokes structure actually occurs at a value of E slightly larger than E" . Hence the quantisation condition associated with this 
arrangement does not come into effect exactly at E = E" 



Table 2: Arrangements of Turning Points, a < 
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(Figure O 
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(Figure O 
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Region 
(Figure [6]l 




E < -E' 


E = -E' 


-E' < E <0 


E = 


< E < E' 


E = E' 


E> E' 


Boundary of E / F 


I = 

















i-E' ^0 for 1 = 1' 
»E' = for I = I' 

* Change in Stokes structure actually occurs at a value of E slightly larger than E" . Hence the quantisation condition associated with this 
arrangement does not come into effect exactly at -E = E" 



To summarise the organisation of these different structures: for a given a, the sequence of 
turning point configurations obtained will depend on the value of / relative to these key values. 
This enables us to divide the (a, /) plane up into sectors according to which sequence of turning 
point arrangements is valid there. A sequence is obtained as the pattern of zeros obtained also 
depends on energy. 

The regions of the {a, I) plane where different turning point configurations occur are shown in 
Figure O This shows how the plane is divided in terms of these several key values of angular 
momentum, /. The boundaries between these regions (shown as solid lines in Figure [H]) are given 
by the lines a = 0, / = and / = ±/'. It is also apparent that the problem ([1]) is symmetric 
about / = — |. 

Note that combinations of these sections approximate sections A - D of Figure [2], and the borders 
of the latter are included in Figure [6] for comparison (dashed lines). Later, we will compare our 
results with those found by Dorey et al. [1,2], and refer back to Figure El (The reader may also 
note that although these regions do not match up perfectly with those found by Dorey et al., 
all results are consistent; using our method, reality and positivity are demonstrated for smaller 
regions than in [1,2], so we have a weaker condition overall - this point will be returned to in 
sections [231 and [231) • 




Figure 6: Regions of the a - I plane . The solid lines show how the plane 
is split up in terms of different Stokes structures occurring. The dashed lines 
indicate the boundaries from [1, 2]. 
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Tables ([T]) and (j2j) are organised so that one first looks at the sign of a to decide which table 
to use (Tabled] is for positive values of a), then which of the five rows to use depending on the 
value of I. Each row shows how the arrangement of zeros changes as E is varied. 
Looking along one of these rows then, it can be seen that change occurs at i? = 0, or one of 
four different critical energy values (labelled i?', E'" and Eq). E' and £'0 refer to energies at 
which turning points come together to form the double zeros described earlier, while E" and E'" 
are energy values at which turning points line up. Expressions for each of these critical values 
(as a function of a. and /) have been obtained (though they are mostly unenlightening) , and 
these energies can easily be calculated exactly for any particular values of the other parameters. 

2.3.1 Example of finding WKB quantisation condition 

Figure [7] shows the structure of Stokes and anti-Stokes lines in the case a = 3, / = 0.5 and 
E = 1. I in this case is less than the critical value /' (/' ~ 1.081 for a = 3), so this refers to 
section / in Figure [6] and Table [1] 



In sector S-i, the WKB wave function proportional to 

(F(x))-^e*''(^i'^) 

is sub dominant, where 



u{a, h) 



no no 

I {p{t))-2dt= / {E-v{t)y2dt 

a J a 



and the Xi are the turning points shown in Figure [71 (In the following, dominance or subdomi- 
nance of a solution is denoted by (d) or (s)). We now continue this solution in an anticlockwise 
direction, and take account of new contributions appearing at Stokes lines and the dominance 
changes at anti-Stokes lines: 

(F(x))-ie*^(^i'^) 

Going from (a) to (b) now involves crossing an anti-Stokes line, so the subdominant solution 
now becomes dominant. We will suppress the {P{x))~i term in the remainder of the discussion 
for clarity. 



(b) 



gja;(xi ,x) 



To pass from region (b) to (c), the Stokes line coming from xi is crossed. The dominant term 
remains unchanged, and the subdominant term gains an extra contribution equal to the domi- 
nant coefficient multiplied by the Stokes multiplier for this line (i). 

^iuj{xi,x) _|_ j^^-iu){xi,x) 

(d) (s) 
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This solution is re- written in terms of the WKB solution from turning point X2- 

Qiijj{xi,X2)QiLLj{x2,x) _j_ j^Q-iLj(xi,X2)Q~ii^{x2,x) 

(d) (s) 

Next a Stokes line from X2 is crossed, and so the subdominant term again picks up a contribution 
proportional to the dominant term. 

QiiLi(xi,X2)Qiuj{x2,x) _j_ ^ ^Q-iuj{xi,X2) _j_ gia;(a:i,a;2)l Q-iuj{x2,x) 

'''^ (d) (s) 

This is again written in terms of the WKB solution from X3, in preparation for crossing the 
Stokes line from that turning point. 

Qiui{xi,X2)Qi'^{x2,x-})Qiui{x-j,x) _|_ i^Q-ii^(xj ,X2) ^Qii^(xi,X2)'^Q-i'^(x2,x:^) Q-ii^(x^,x) 

^ (d) (s) 
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This procedure repeats, as a further Stoke hne is crossed. 



(d) (s) 
(d) 

_|_ j| |^e-»"(a:i,a:2)_|_g»"(a:i,a:2)jg-»"(a;2,a;3)_|_giij(a:i,a:2)e'"(^2i^3) j.e~'"(^3'^4)e~'"(^4i^) 

(s) 

/r\ Qi^(xi.X2)Qiui(x2,x^)Qiui(x3,X4)Qiui(x4^,x) 

(d) 

[G~''^'^i'^2)+e''^f^i'^2)je~*"(^2.^3)+e''^(^i'^2)e''^(^2.^3) j.e-«'^(^3>^4) 
_i_ e»"(^i.^2)e'"(^2.^3)e'"(^3'^4)je~'"^'^4'^) 

Finally, in going from (f) to (g), an anti-Stokes line is crossed again, and dominance is again 
exchanged. 

Qiu}(xi,X2) QiLij(x2,X3) Qiu!(x3,X4,) Qiu!{x4,x) 

(s) 

_|_^ [{ l^e"*'^^^^'^^-' + g*'^(^l'^2)j g-jt^(a:2,a;3) _j_ gja;(xi,a:2)giaj(x2,X3) | Q-iuj{x3,X4) 
_j_ giaj(xi,X2)giaj(x2,X3)gjLj(a:3,X4) j g-ia;(a;4,x) 

(d) 

So for the solution to be purely subdominant at (g), (Sector Si), in order to satisfy the boundary 
condition, we require 

I ^^-iLj{xi,X2) _|_ ^iiLi{xi,X2)~^ g-ia;(x2,X3) _|_ gjtj(a:i,X2)gia;{a;2,a:3) | ^-iLLi{x3,X4) _|_ ^iiLi{xi,X2)Qiuj{x2,X3)^iuj{x3,X4) _ q 



Writing uj{xi, X2) = U + iV = (a;(x3, X4))* and u{x2, X3) = W, with ?7, F and W real, (a;(x2, X3) 
can be seen to be purely real as X2 and X3 are joined by an anti-Stokes line), following the 
example of [17], this simplifies to: 

E eigenvalue 2 cos(2;7 + W) + 2e-^^ cos{W) = 

As explained earlier, the form of the quantisation condition changes, depending on the parame- 
ters a and /, as well as changing as energy is increased moving along a sequence. We now define 
(piecewise), an overall quantisation condition function, which we will denote Q{E,a,l). e.g. as 
we have just seen Q{E,aJ) = 2 cos(2[/ + VI/) +2e-2^ cos(Vr) for < < , a > , < / < /'. 
Similarly Q{E, a, I) is defined to be the quantisation condition obtained for each of the other 
different sets of parameter values. Finding eigenvalues is then a matter of finding the zeros of 
this function. 
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2.4 Positivity of Spectra from WKB Quantisation Conditions 



The sections E and F in Figure [6] form an area that is completely contained in the sector where 
positivity was proved in [1] (bounded by the red lines in Figure [6l or section D in Figure [2]). 
WKB analysis of E and F shows that no quantisation condition can exist for negative energies. 
In fact, a viable quantisation condition is only found for E > E'. For energy values lower than 
this, demanding that the dominant wave function vanishes in the Stokes sector Si also means 
that the subdominant wave function will vanish; there is no way to have a WKB solution that 
is purely subdominant in sectors iS_i and Si, and hence there are no eigenvalues. 

2.5 Conjecture on Reality of Spectra from WKB Quantisation Con- 
ditions 

In [17], the appearance of complex eigenvalues was explained in terms of a WKB quantisation 
condition. In the paper. Bender et al. obtained a WKB quantisation condition for the VT 
symmetric potential 

V{x) = x'^ + iAx. 

The condition found was: 

E eigenvalue <^=^ cos(2t/) + ^e"^^ = 0, 

where U = (^J^ \e - V{t))U?j , V = ^m(^j^ \e - V{t))Ut 

This condition, consisting of the sum of an exponential term and an oscillatory term can lead 
to real eigenvalues, provided that the magnitude of the exponential term does not exceed 1. If 
parameters are varied, leading to an increase in magnitude of the exponential term beyond 1, 
then no real solutions are possible. In effect, real eigenvalues change in value as the parameter 
A is varied continuously, then 'pair off' and disappear as the exponential term becomes larger. 
This can be seen explicitly by plotting the value of the quantisation condition function against 
energy for a particular value of parameter A, then looking at how this plot changes with different 
values of A. As the value of the exponential term increases with changing A, local minima in the 
plot can be seen to be pulled up through zero, leading to degenerate eigenvalues and complex 
conjugate eigenvalues, (see Figure [S]). 

The WKB quantisation conditions found in regions / and J of Figure O (which are slightly 
larger than the region A (Figure [2]) discussed in [1,2]), were found to be similar in form to that 
obtained in [17] {i.e. containing an exponential and oscillatory term). 

With these conditions the effect described above can occur, where pairs of real eigenvalues move 
together and become degenerate as the parameters change, due to the effect of the exponential 
term(s) in the quantisation condition. 

An example of this is seen by looking at the quantisation conditions obtained for Region / (Table 
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(a) 





Figure 8: Quantisation condition function plotted against energy 



8(a) 



corresponds to 2 real eigenvalues, 8(h) to real, degenerate eigenvalues and 8(c) 
to complex conjugate eigenvalues. (Referring to two lowest energy eigenvalues 
shown). 



[T]). In section [2.3.11 we showed that the quantisation condition for region / (a > 0, < / < /') 
with < E < E' is 

E eigenvalue 2 cos(2?7 + l^) + 2e-2^cos(iy) = 

Plotting this quantisation condition function against energy, for various values of a and /, it 
can be seen that varying the parameters can again lead to the exponential term increasing in 
magnitude so that it begins to dominate the oscillatory term, and pulls the function through 
zero. One important difference to [17] is that the exponential term in the condition is multiplied 
by a cosine term. This makes it harder to analytically pin down the boundary in the parameter 
space where degenerate eigenvalues may occur, than in [17]. (Hence we can only restrict the 
domain of unreality to the regions where particular quantisation conditions are valid, and don't, 
as yet, make a more detailed study of these regions - i.e. we only say that complex conjugate 
eigenvalues exist somewhere in this region, and do not analytically calculate the particular area 
inside this region in which real eigenvalues are not possible (c.f. [17])). 

This cosine term leads to another difference in the quantisation condition function plots; because 
this term varies between ibl, the exponential term in the condition will be positive for some val- 
ues of a, I, E and negative for others, degenerate eigenvalues can now come about through a 
local minimum being pulled up through zero, or a local maximum being pulled down. See Figure 

IS 



In all other regions, the quantisation conditions do not contain exponential terms. The most 
common quantisation condition found outside of I and J is of the form cos([/) = 0, apparently 
leading to an infinite number of real eigenvalues. We conjecture that the only complex eigen- 
values found for this problem are the ones that come about in the manner described earlier, 
whereby the exponential term allows for 'pulling through' of the quantisation condition func- 
tion. If it is true, as seems to be the case numerically, that this is the only means of complex 
eigenvalues coming into existence, then any quantisation condition lacking an exponential term 
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Figure 9: Quantisation condition function, illustrating degeneracies arising 
from mimimum being pulled up, and maximum being pulled down. 



(or at least a sum involving two similar sized terms) will not lead to complex eigenvalues. Hence 
our conjecture leads to the conclusion that, in Figure E], the spectrum will be entirely real below 
the blue lines in the upper part of the plane. This is a weaker condition than that found in [1], 
where reality was established for all parameter values below the red lines in this upper part of 
the plane. 

(Note that there is one exception to the statement that the quantisation condition does not 
contain exponential terms outside of I and J; that being the case I > 0, E > E", i.e. high 
energy values in regions E, G, H and /; The quantisation condition found here is of the same 
form, 2 cos(2[/ + W) + 2e~^^ cos(VF) = 0, but this condition only comes into effect for E > E" 
whence it is seen numerically that V is relatively large and positive, so the exponential term 
is approximately zero and the pairing off of eigenvalues cannot occur, and so again only real 
eigenvalues are possible). 

2.6 WKB explanation for appearance of degeneracies and cusps 

The pattern of cusps in the degeneracy curves described in [2] is manifest from the WKB ap- 
proach via the quantisation conditions. The position of degenerate eigenvalues can be found 
numerically by searching for the parameter values where a local maximum (minimum) in the 
quantisation condition passes through 0. this is done by combining a simple maximum (mini- 
mum) finding algorithm, and a bisection algorithm to find when the maximum (minimum) value 
of the quantisation condition function is 0. In this way, a plot of the locations of degenerate 
eigenvalues can be found. (Figure [TOj) . 

Figure dU] shows the pattern of curves, meeting at cusps, as seen in Figure [31 first shown by 
Dorey et al. [1]. Note that this figure is different to Figure [31 as it shows the cusp pattern 
repeating itself as further eigenvalues join up and then split into complex conjugate pairs. The 
results from the complex WKB quantisation conditions are shown as the crosses, overlaying the 
results from the method used in [1,2]. It can be seen that there is extremely good agreement 
between the degeneracy curves found via these two methods. We are extremely grateful to 
Patrick Dorey, Clare Dunning, Anna Lishman and Roberto Tateo for sharing their results prior 
to pubhcation [18,19]. 

To explain how the quantisation condition method accounts for the formation of the cusps in 



16 



-10 -5 5 10 

/ 

Figure 10: Degenerate eigenvalues The positions of degenerate eigenvalues 
found using WKB quantisation conditions (crosses) and those found with the 
methods of [1,2] (lines). 



the diagram, consider the diamond-shaped area of the (a, /) plane bordered by < «„ < 1 
and < < 1, In this region, for a fixed value of «_)., and a_ being increased from to 1, 
degeneracies are seen as a local minimum of the quantisation condition function being pulled up 
through zero. In the adjoining diamond-shaped area, < a_ < 1 and 1 < a_|. < 2, degeneracies 
occur as a local maximum is pulled down through zero. In < «_ < 1 and 2 < a+ < 3, it is 
again a minimum, and so on (Figure [TT]). 

On the boundaries of these regions, ie a± G Z"*", the degeneracies occur as a limit of a local mini- 
mum and local maximum of the quantisation condition function both meeting at Q{E., ct, /) = 0, 
ie a point of inflection of the quantisation condition function. These correspond to the cusps 
in [2]. 



An example is shown in Figure [T21 12(a) shows the behaviour of Q with parameter values ap- 



proaching those for which a point of inflection occurs. It can be seen that both a local maximum 
and local minimum are close to Q = 0. Figure 12(b) corresponds to a degeneracy for slightly 
less than 2, (a_ ~ 1). Here the minimum has just been pulled up to Q = 0. 12(c) corresponds 
to a degeneracy for a+ just greater than 2, («_ ~ 1 still). Here the maximum has just been 
pulled down. 12(d) shows Q for the parameter values a+ 2,a_ 1, leading to an inflection 
point. 

Numerically, by scanning the (a, /) plane for points of inflection, a plot of the positions of these 
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Figure 11: Degeneracies being formed from minima or maxima of 

QiE^ a, I) On the edge of one diamond shaped region, degeneracies occur as 
either a minimum or maximum. In the neighbouring diamond this is swapped. 




(a) (b) (c) (d) 

Figure 12: Quantisation condition function, parameter values close to 
inflection point. 



cusps, as predicted by WKB was obtained (Figure [T3l). 

These cusp positions are in approximate agreement with those found by Dorey et al. in [1,2]. 

3 Conclusions and Further Directions 

We have shown that the approach of complex WKB quantisation conditions can be used to 
explain aspects of the spectrum of ([1]). The problem appears to be more complicated due to 
the changing of the Stokes structures as parameters are varied, but in fact it is this changing 
of configurations that provides an explanation for the different behaviour in different regions of 
the (a, /) plane. 

The M = 3 case is a useful test case to investigate as some energy levels can be found exactly, 
due to the quasi-exactly solvability of the related problem discussed in [2]. The quantisation 
condition method has been shown to give excellent agreement in calculating the positions of the 
degenerate eigenvalues, (possibly surprisingly good, as first order complex WKB is supposed to 
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Figure 13: Locations of cusps. 



be a more accurate approximation for higher energy levels, and the pairing off occurs between 
low lying energy levels). There should be no reason why the method would not be useful 
in investigating the spectrum for different values of M. Indeed, considering smaller values of 
M would be expected to be a simpler problem, as less zeros of the potential should lead to 
less complexity in the pattern of Stokes structures obtained. It may be possible to see how 
the quantisation condition approach changes for general M, and gain insight into the way the 
domain of unreality changes as M is varied. 
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